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We provide a detailed quantitative analysis of transport properties in the p-wave supercon- 
ducting state of Sr2RuC>4. Specifically, we calculate ultrasound attenuation rate and electronic 
thermal conductivity within the mean field approximation. The impurity scattering of the 
quasi-particles are treated within the self-consistent T-matrix approximation, and assumed to 
be in the unitarity limit. The momentum dependence of the gap function is determined by 
solving the Eliashberg equation for a three-band Hubbard model with the realistic electronic 
structure of S^RuCU. On the basis of the microscopic theory, we can naturally expect nodal 
structures along the c-axis on the cylindrical Fermi surfaces, even if we assume the chiral pair- 
ing state (i.e., A(k) ~ k x ± ifcy). Consequently, we obtain the temperature dependence of the 
transport coefficients in agreement with the experimental results. We can clarify that actually 
the thermal excitations on the passively superconducting bands contribute significantly to the 
thermal conductivity in a wide temperature range, in contrast to the case of other physical 
quantities. 
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trasound attenuation rate, thermal conductivity 



1. Introduction 

The quasi-two-dimensional ruthenate Sr2Ru04 has at- 
tracted much interest of solid state physicists since the 
discovery of superconductivity in this compound. 1 A 
number of excellent experiments and theoretical consid- 
erations have already revealed its remarkable and in- 
triguing physical properties. The most important is that 
Sr2RuC>4 is a strong candidate of odd-parity spin-triplet 
superconductor (most likely p-wave superconductor), 1-3 
although the crystal structure of Sr2Ru04 is the quasi- 
two-dimensional pcrovskitc structure identical to that of 
high-T c copper oxides (spin-singlet d x 2_ y 2-wetve super- 
conductors). 

Mechanism of the spin-triplet pairing in Si'2Ru04 
has been discussed from various microscopic theoretical 
points of view. In the early stage of the researches, it was 
considered that some strong ferromagnetic spin fluctua- 
tion (Paramagnon) will exist and induce the triplet pair- 
ing. ■ 5 However, the inelastic neutron scattering mea- 
surement showed that the spin correlation is predom- 
inantly incommensurate antiferromagnetic rather than 
ferromagnetic. 6 Interestingly, some recent theoretical 
works have shown that the pair-scattering amplitude 
could actually possess such a characteristic momentum 
dependence as is favorable for p-wave pairing but not 
attributable to that of the spin susceptibility. 7-9 These 
theoretical discussions suggest that the triplet pairing in 
Sr2RuC>4 is considered as a natural result from electron 
correlation, but the essential pairing attraction between 
the Fermi liquid quasi-particles does not originate from 
the exchange processes of spin fluctuations. 7-9 

Another interesting problem, which is attacked in the 
present work, is the detailed clarification of the supercon- 
ducting gap structure. Superconducting gap structure is 
closely related to the temperature dependence of phys- 



ical quantities in general, since the thermal excitations 
only near the gap minima are essential in responding to 
various external perturbations. On the theoretical side, 
the gap structure is a reflection of the momentum depen- 
dence of pairing interaction. Therefore, the detailed anal- 
yses of the temperature dependence of physical quanti- 
ties could be a valuable test for the validity of pairing 
scenario. 

The gap structure of Sr2Ru04 is still controversial. 
The most probable pairing symmetry A(fc) ~ (k x ±ik y )z 
was proposed by Rice and Sigrist very early after the 
discovery of the superconductivity. 4 The muon spin re- 
laxation experiment 10 suggests strongly that this chiral 
pairing state, in which the time reversal symmetry is 
broken, is the most plausible pairing state among sev- 
eral candidates. This pairing symmetry usually results in 
isotropic gap (i.e., the absence of nodal structure). How- 
ever, most of existing experimental data exhibit power- 
law temperature dependence, 11-16 and indicate a nodal 
gap structure (most probably line nodes) on the Fermi 
surface. Much effort has been devoted to reconciling these 
contradictions. 17-23 

In the present article, we analyze microscopically the 
temperature dependence of two transport coefficients, ul- 
trasound attenuation rate and thermal conductivity, in 
the superconducting state of Sr 2 Ru04. The sound atten- 
uation rate is expected to provide a strong constraint on 
the theoretical proposals of superconducting gap struc- 
ture. 24 We derive the momentum and band dependence 
of the gap function by solving Eliashberg equation for 
a realistic tight-binding electronic structure. 23 Starting 
with a three-band Hubbard model, the effective pairing 
interaction is evaluated perturbatively to third order in 
the on-site Coulomb integrals. 8 In our previous works, we 
showed naturally by this theoretical framework that (i) 
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the momentum dependence of the pairing interaction is 
favorable for anisotropic p-wave pairing, and p-wave state 
is obtained as the most probable pairing state, 7 ' 8 (ii) 
the superconducting gap possesses line-node-like struc- 
tures along the c-axis, 9 ' 23 and (iii) the calculated specific 
heat as a function of temperature fits experimental data 
well. 23 

The transport coefficients of Sr 2 Ru0 4 have been cal- 
culated by other authors. 25-30 Most of them are based 
on simplified isotropic Fermi surfaces (two- or three- 
dimensional), or simple gap functions described by one 
or a few harmonic functions. However, as discussed in 
the present article, such simplified electronic structures 
and gap structures are generally insufficient to discuss 
the experimental data of transport properties. 

For the calculation of the transport coefficients, we 
take the way of analysis which has been developed 
for some uranium compound superconductors: 31-33 the 
non-magnetic impurity scattering is treated by the self- 
consistent T-matrix approximation. The scattering is as- 
sumed to be in the unitarity limit. We ignore the ver- 
tex corrections. In contrast to the analyses for uranium 
compounds (where they considered simplified electronic 
structures, e.g., isotropic spherical Fermi surfaces), wc 
take more realistic tight-binding electronic structure and 
gap structure for S^RuO^ 

For the calculation of ultrasound attenuation rate, we 
adopt the same way as Walker and collaborators did. 34 
In the work by Walker et al, they succeeded in explain- 
ing the large in-plane anisotropy of the attenuation rate, 
and clarified that the anisotropy of the electron-phonon 
interaction is essential. This successful point is taken into 
account also in our present work. However, it should be 
noted that their gap structure is quite different from 
ours. They insisted that there should be point nodes on 
the Fermi surfaces. 30 Although we discuss based on the 
purely two-dimensional model, we consider that it is al- 
most impossible to obtain point nodes or horizontal (i.e., 
parallel to the basal plane) line nodes. This is because 
the pair scattering amplitude can acquire only negligi- 
ble momentum dependence along the c-axis, due to the 
strong two-dimensionality. 

Consequently, we obtain a consistency with the exper- 
imental data of the transport coefficients in the overall 
temperature range. The electron-phonon coupling con- 
stants are essential parameters for the good fitting to the 
experimental data. In addition, the anisotropy due to the 
lattice structure is also an essential ingredient for repro- 
ducing the thermal conductivity and the anisotropy of 
sound attenuation. Therefore the naive conjectures from 
the isotropic models and simplified gap functions are not 
reliable in general. 

The present article is constructed as follows. In § 2, 
the theoretical formulation is given. We give the analytic 
expressions of ultrasound attenuation rate and thermal 
conductivity. In § 3, the numerical result of calculated 
gap structure is provided, and then the numerical results 
of the ultrasound attenuation rate and thermal conduc- 
tivity are compared with the experimental data. In § 4, 
some remarks on our results are given. In addition, rela- 
tions of the present work with other works are discussed. 
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Finally, in § 5, we will give concluding remarks. 
2. Theoretical Formulation 

2.1 Three-band Hubbard model and Eliashberg equation 
The Fermi surface and the electronic structure of 
Sr2Ru04 near the Fermi level are well reproduced by 
tight-binding fitting. 35 Since the electronic density of 
states near the Fermi level is dominated by the partial 
density of states of Ru4ds orbitals, 35 ' 36 the hopping pa- 
rameters are considered to describe the transfers between 
the Ru4de-like Wannier orbitals. We take the following 
non-interacting Hamiltonian (these Wannier orbitals are 
characterized by I = {xy, yz, xz}): 

+ A ( fe ) ( C fe If* o C k xza+c{ xza C kyza ),(l) 
h.cr 

where Ckia[c kicr } is the electron annihilation[creation] op- 
erator (the pseudo-momentum, orbital and spin states 
are denoted by k, I and a, respectively), and the energy 
dispersions are 





= (cos k x + cosfc y ) 






+4^2 cos k x cos k y — n xy , 


(2) 




= 2^3 cos ky + 2^4 cos k x — fj, yz , 


(3) 




= 2t 3 cos k x + 2t 4 cos k y — fi xz , 


(4) 


A(fc) 


= sin k x sin k y . 


(5) 



In the present work, we take the parameter set t\ = 
-1.00, t 2 = -0.395, t 3 = -1.25, U = -0.125, t 5 = 
—0.150, 37 to reproduce the Fermi surface topology. The 
chemical potentials fig's are determined by the condition 
ni = 0.667, where ne is the electron number of orbital 
£ per one spin state (The orbitals arc evenly filled with 
electron). We introduce the Coulomb interaction part: 

+y EEE 4<t4'<t' c ^vq£ CT 

i t=£l' ct,<t' 

+ yEEE cl^cl^CWa'CWa, (6) 

where the operator cu a [c\ t(J ] is the electron annihila- 
tion [creation] operator at i-th Ru site (ckia is the Fourier 
transform of cu a ). The microscopic origin of H' is the 
Coulomb interaction between the Ru4d electrons. The 
total Hamiltonian is H = Hq + H' . 

As we will see in the next section, the Fermi surface 
consists of three branches (which are named a, (3 and 7), 
consistent with the de Haas- van Alphen oscillation 38 and 
photoemission measurements. 39 The Hamiltonian H is 
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diagonalized easily, and the obtained dispersions are 

Uk) = z + (k)-y/e_(k) + x>(k), (?) 

fr(fc) = Z + (k) + y/£(k) + X>(k), (8) 

£ 7 (fc) = &„(fc), (9) 

where £±(fc) = (£y Z (k) ± S, xz (k))/2. The elements of di- 
agonalization matrix arc 



U(k) 





1 


a 


P 


xy 


1 








yz 





K{k) 


L(k) 


xz 





-L(k) 


K{k) 



(10) 



where K{k) = ^(1 - M(k))/2, 

L(k) = sgn(A(fe))V(l + M(fe))/2, and M(fe) = 

£-(*0/\A-(*0 + A 2 (fc). We define the bare Green's 
function: 



G^(k) = (ico n -Uk)Y 



(11) 



the band index a takes a, (3 and 7. 

The anomalous self-energy (i.e., superconducting order 
parameter) is obtained by solving the linearized Eliash- 
berg equation: 



T 



a' ,k' ,C3CT4 



(k,k') 



x\GZ>(k'WK>,w 3 (n (12) 

where fcW = (fc (,) , iw„ (/) ) [w n = (2n+ 1)ttT], and band 
index takes a, /3 and 7. The effective interaction 
Vaa 1 a 2 .a'a 3 cr 4 {k, k') is evaluated by the third order per- 
turbation expansion in H'. In the present calculation, 
we take U = 4.0, U' = 0.33L7, J = 0.33L7, J' = 0.3317. 
The lengthy procedure of expansions and the numerical 
solutions of the Eliashberg equation (12) are given in 
the previous work. 8 For the spin-triplet states, the or- 
der parameter is expressed using the vectorial notation: 
^, aia2 (k) = [i(D a (k) ■ <j)a y ] aia2 . 

We assume that the momentum and temperature de- 
pendence of the superconducting gap function A a (fc) arc 
given by the following equation (as in ref. 23), 

A„(fc) = f a (k)A(T), (13) 

where f a (k) = Const, x D a (k), and D a (k) is the compo- 
nent of the vector D a {k) (we assume D a (k) = D a (k)z.), 
and obtained by solving the above Eliashberg equa- 
tion (12). The temperature dependence of the gap mag- 
nitude A(T) is determined by the standard BCS gap 
equation: 

tanh[£ a ,(fc')/(2T)] 



a'k' 



E a ,{k') 



With Vfca.fc'a' 



L A ,(fc'), 
(14) 



-f a {k)f*,{k') and E a (k) 



V£a(k) + |A a (fc)| 2 - fa(k) is an odd-parity function, sat- 
isfying the relation f a (k) = —f a (—k). We assume the 



chiral state is realized: 

/a(fc) = / *(fc)±i/»(fc), 



(15) 



where fa{k){f%(k)] is a real function possessing the k x - 
like[fc y -like] symmetry. The superconducting gap struc- 
ture on band a is obtained by the absolute magnitude of 
A„(fc), 

|A„(fe)| = |/„(fc)|A(T) = y/fZ(k) 2 + f!>(k) 2 A(T). 

(16) 

For the details of the discussions in this section, one 
could refer to refs. 8 and 23. 

2.2 Self- consistent T -matrix approximation for non- 
magnetic impurity scattering 
Here we consider non-magnetic impurity scattering, as 
preliminaries for the following discussions on transport 
properties: 



H-, 



imp = Y ^ J U it lC jla C 3t'<y- 
J Gimp. W a 



(17) 



"j G imp." means that the summation is performed over 
the impurity positions. The impurities are randomly po- 
sitioned at Ru sites. We assume simply uw — u5a>, be- 
cause of the symmetry properties of the Ru4de orbital 
wave functions. In addition, we assume that u is much 
larger than the band width W (u > 100W for the nu- 
merical calculations), i.e., the scattering is almost in the 
unitarity limit. 

In order to discuss the effect of impurities, we take the 
standard T-matrix approximation. 40 The self-energy for 
the impurity scattering is given by 

£«/(iw„) = cT u ,(iw n ). (18) 

c is the impurity concentration, and we take c = 10~ 6 
for the present work. Tw(iu) n ) is the T-matrix obtained 
by solving the equation, 

T«/(iw„) = u U ' + -jy Y u «iG^ 2 (fci,iw„)T£ 2 Wiw n ), 

fel,«2 

(19) 

where Gw(k) is the Green's function. We have already 
omitted the terms containing X^fc-^'W [Fw{k) is the 
anomalous Green's function. See below.], which vanish in 
non-s-wave superconducting states, 42 in contrast to the 
case of s-wave superconducting state. The impurity self- 
energy and Green's functions are related by the Gorkov 
equation: 

G a {k) = G(°)(fc) + G£ >(fc)£ o (fc)G o (fc) 

-Gi°)(fc)A a (fc)F a t(fc), (20) 

F a (k) = G<°>(fc)£ a (fc)F a (fc) 

+Gi°)(fc)A a (fc)G a (-fc), (21) 

Fl(k) = Gi°)(-fc)E a (-fc)F a t(fc) 

+G^(-k)A* a (k)G a (k). (22) 
Using the matrix U(k) in cq. (10), the orbital indices &^ 
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in E«'(iw„) are converted to band index a by 

E (fc) = ^C/^(fe)^ (fe)S«,(io;„), (23) 



and the band index a in the Green's functions are con- 
verted to the orbital indices &^ by 

Gu>{k) = ^C/, Q (fc);7 a t ,,(fc)G a (fc), (24) 

a 

Fw(k) = Y,Uta(k)U ea (k)F a (k), (25) 

a 

iMfc) = E^W^'W^W- (26) 

a 

Defining the rcnormalized Green's function by 
G a (k) = (Gi )" 1 ^)-^))- 1 

= (iWn-^fcJ-E^fc))- 1 , (27) 
the above Gorkov's equation is simplified to 

G a (k) = G a (k)-G a (k)A a (k)F^k), (28) 
F a (k) = G a (k)A a (k)G a (-k), (29) 
Ft(fe) = G„(-fe)A:(fc)G„(fc). (30) 

By solving eqs. (28)-(30), the Green's functions are ob- 
tained: 

1 



G a (k) 
F a (k) 
Fl(k) 



GaW + GKWAaiW 

Aq(fc) 

|G a (fc)|-2 + |A a (fc)P' 
A*(fe) 



(31) 
(32) 
(33) 



|G Q (fc)|- 2 + |A a (fc)| 2 ' 

We continue this expression analytically to the real fre- 
quency axis (ko n — > z + iO): 

G R (fc) = 

^(-fc)-ea(fc) 

(~z*(k) - Uk))(^*(-k) &(*)) + |A a (fe)| 2 ' 

(34) 



F*(k) 



A a (fe) 



(zTW - e a (fc))(^*(-fc) - e a (fe)) + |A a (fc)i 2 ' 

(35) 

^a tR (fc) = 



(z*(k) - Uk))(^*(-k) Uk)) + |A a (fe)| 2 ' 

(36) 



with 



(37) 



and = (fc,z). Here we may assume the electron- 
hole symmetry in order to simplify the discussion, al- 
though the actual electronic structure does not possess 



the symmetry. This simplification does not affect the re- 
sults, since the temperatures we consider here are enough 
low, compared with the characteristic energy scale of 
the asymmetric structure of the density of states. Un- 
der the electron-hole symmetry, we have the relation 



to 



(-*0 



zR 



„ (k), and the expressions are simplified 



G R (fc) 



s«(k) + Uk) 



~z*{kf-Uk? 



|A a (fc)P 



F*{k) = — — 



A a (fc) 



z^{kf-Ukf 



|A a (fe)| s 



Fl R (k) = - 



K(k) 



(38) 



(39) 



(40) 



S^(k) 2 -Uk) 2 -\A a (kW 

We determine z n (k), by solving simultaneously the an- 
alytically continued forms of eqs. (18), (19), (23) and 
(24), and eqs. (37), (38). There we use A Q (fc) derived 
in the absence of the impurities (in § 2.1). This means 
that the changes of the gap function and its temperature 
dependence due to the impurities are neglected in the 
present study. We consider that the impurity concentra- 
tion c(= 10~ 6 ) is enough small not to change significantly 
the gap structure and its temperature dependence. 

In the present work, we assume that the damping of 
thermally excited quasi-particles is dominantly caused by 
the impurity scattering, and ignore the damping effect 
due to electron-electron scattering. We could justify the 
assumption as follows. Within a simple discussion, we 
expect the normal-state electronic thermal conductivity 
n n behaves at low temperatures as 

n- 1 ~ dT- 1 + fT + gT 2 , (41) 

where T is temperature, and d, / and g are constants. 41 
The first, second and third terms of the right hand 
side in eq. (41) originate from the impurity scattering, 
the electron-electron scattering and the electron-phonon 
scattering, respectively. According to the experimental 
data of the normal-state thermal conductivity in ref. 16, 
we find that the first term is sufficiently larger than the 
other two terms within the low-temperature region. 

2.3 Ultrasound attenuation rate 

Scattering of phonons by the electron system causes 
the ultrasound attenuation. We consider the electron- 
phonon interaction: 

H ep = N~5 E ^k,q,U"t>qCk +q e< 7 Ckt>a, (42) 
kq,U',a 

where 4> q = b q + b^_ q , and b q (b q ) is the phonon annihila- 
tion(creation) operator with momentum q. The matrix 
elements Afe i9i «/ for the present electronic structure of 
Sr 2 Ru04 are given by 

A-k,q,xy,xy ~ '^\9l (COS k x C x q x ~t~ COS kyCyQy ) 

+g 2 cos k x cos k y (e x q x + e y q y ) 
-g 2 smk x smk y (e x q y + e y q x )}, (43) 
•U.,,.., = i{g4 cos k x e x q x +g 3 cos kye v q y ), (44) 
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G(k) 



G(k+q) 



-k 



F(k) 

-k-q 3 V v> ____ < ^>k+q 

F f (k+q) 



Fig. 1. The diagrammatic representation of the phonon self- 
energy part within the mean field approximation. The solid lines 
with arrowheads denote the Green's functions. 



Ak,q, X z,xz 
Ak,q,yz,xz 



i(<7 3 COS k + gA cos k y e y q y ), (45) 

Ak,q,xz,yz 

= 155 {- sin fez smk y (e x q x + e y q y ) 
+ cos k x cos k y (e x q y + e y q x 

)}, (46) 

and the other elements are zero (See Appendix A), 
e = (e x ,e v ) and q — (q x ,q y ) are the unit vectors along 
the directions of phonon polarization and phonon prop- 
agation, respectively. 

The energy spectrum of phonon is given by the poles of 
phonon Green's function. The phonon Green's function 
D{q) [q = (q, if^m)) = 2mirT] satisfies 

D(q) = D(°\q) + D^(q)n(q)D(q), (47) 

where D^°\q) is the bare phonon Green's function 
[£><°>(«) = -2cu (q)/(n m 2 +uj Q (q) 2 )}, and U{q) is the 
phonon self-energy. Using the retarded phonon Green's 
function D R (q) [D R (q) is the analytic continuation of 
D(q) to the real axis by ifi m — > Q + iO (m > 0)], the 
attenuation rate a of the phonon with momentum q is 
obtained by solving D R (q,Q) = for ft and then 
taking the imaginary part of the solution: a = Imfl f» 
— lmIl R (q, ojo(q)). H R (q, f2) is the analytic continuation 
of n(q, ifl m ). The phonon self-energy H(q) within the 
mean field theory is represented by the diagrams in 
Fig. 1. The analytic expression of the phonon self-energy 
is 



Mi 



Ak.q.l 



q,C 3 U 



x[Gt 2 £ 3 (k + q)G tiil (k) 



q)F tA {k% (48) 



where we have used the relation A_k, q ,U' = ^-k,q,W, 
and the approximate relation Ak+is q , q ,u' ~ Afe.q,«' (for 
v ~ 1) within the precision to the leading order of |q|. 
The front factor two originates from the summation with 
respect to spin indices. We continue II (q) analytically to 
the real axis, convert orbital indices ii to band index a, 
and take the imaginary part. Since we consider the hy- 
drodynamic limit (£ c <C 2n/\q\, £ c : mean free path) in 
the present study, we may take the limit q, £1 — ► in the 
argument of the Green's functions. Thus we obtain the 
expression for the attenuation rate a: 



^o(g) 

2ttT 



/oo 
dz 



1 



cosh 2 (^) 



x[(ImG*(fc,z)) 2 - |A a (fc)| 2 (ImF a R (fe,z)) 2 



(49) 



where F*(k, z) is defined by Ff(k, z) = A a (k)F^(k, z), 
and wc have neglected the inter-band crossing terms, 
ImG a (k)lmG a ,(k) and ImF a (fc)Im^(fc) (a ^ a'). The 
contribution from these crossing terms is considered to 
be negligibly small in the present case of the long- 
wavelength limit and low temperatures. Ak, q , a is related 
to A k , q ,u< by 

A k ,q,a = u liS k ) u t2a{k)A Kq4ll , 2 . (50) 

lit? 

Here wc perform the momentum integration perpen- 
dicular to the Fermi surfaces in eq. (49). Using the ex- 
pressions of the Green's functions (38) and (39), and 
adopting the approximate form of integration measure, 



£ 

a,k 



£ 



dk 



1 f 

k=k F •> -< 



dCa 



the formula (49) is reduced to 



a 



2ttT 



1 



cosh (2^) 



kp ,q,a\ 



Im- 





— 1 pOG 


dk 


k=k F J -co 



(51) 



d£a 



z a R(fc F ,z) 2 -e a 2 -|A a (fe F )|^ 



Im- 



|A„(fci 



zH(fc F ,z) 2 -C a 2 -|A a (fc F )P 



2n 



(52) 



X)fc F means the momentum integration on the Fermi sur- 
face. We obtain the final expression of a by performing 
the integration with respect to £ a : 



a = 



up{q) 
8T 



with 



dZ u2/ z Ia(k F ,z), 

-00 cosh (2%) 

(53) 



dk 



fe=fe F 



I mv /zR(fc F ,z) 2 -|A a (fe F )| 2 

1+ \z« {kF ,z)f-\A a (k F) n 

|z a R(fc F ,z) 2 -|A a (fc F )l 2 |/ 
The electron-phonon coupling A& i9i(J must be renor- 
malized to satisfy the local charge neutrality condition 
(See Appendix B, and also ref. 41): 

A k 



Ak,q,a Aq, 



(55) 



where A q is the mean value of A/-,q, a over the Fermi 
surfaces, and given by 

Ea, fc {(%G)a(fc)}AM,q 
Ea, fc (^G) a (fc) 

E„ AV2 {(^G)a(fc)}<(fe)t/^(fe)A M , W2 



A q = 



Ea, fe (^G) a (fc) 



(56) 
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and 



Ag(fc) 2 

2E a (kf 



tanh 



\ 2T 



4rK(fc) 2 cosh 2 (^l) 



• (57) 



Throughout the study of sound attenuation, we always 
retain the condition, and use the renormalized h-k,q,a- 

2.4 Thermal conductivity 

The thermal conductivity tensor is calculated employ- 
ing Kubo formula: 42-44 

K "" = ¥ 75 Um n ImP ^(«' )' ( 58 ) 

1 i I q^O ^ 

where P^ v (q,fl) is the retarded thermal-flux correla- 
tion function, and obtained by continuing analytically 
m ) is the Fourier transform of 

P„„(<f,r) =J2< TrJ^(q,T)jl(-q,0) >, (59) 

a 

where r) is the imaginary-time thermal- flux oper- 

ator on band a: 



k,a 

xcLWcfe+^^r'), (60) 
and Wfc p is the ^-component of the velocity on band a: 

= -^r- (61) 

Within the mean field theory, the correlation function 
P(q, i^m) is expanded using the Green's functions as 

P MI /(g,ifi m ) = - — y^{(^n + n m )vkan + ^n«fc+qa M } 
a.fc 

X{(w„ + fl m )Vkav + ^nVk+qav} 

x[G a (k + q)G a (k) - Fl(k + q)F a (k)]. (62) 

We perform the analytic continuation if2 m — > Q + iO of 
cq. (62), and substitute it into eq. (58), then we obtain 
the expression for the thermal conductivity n^ v : 

1 f°° z 2 

K ^ = 2*F 2 ^ J-°° dZ cosh 2 (#) 

x[(ImG^(fc,z)) 2 -|A a (fc)| 2 (ImF a R (fc,z)) 2 ]. (63) 

We adopt the approximate integration measure in 
eq. (51), and perform the integration with respect to £ a , 
as we have done for the sound attenuation rate (in § 2.3). 
We obtain the final expression for the thermal conduc- 
tivity: 

1 f°° z 2 
K ^ = gjv Yl v kFa^k Fav J dz u2/ z J a (k F , z), 



3. Numerical Analyses of Experimental Data 

3. 1 Gap structure and density of states 

In the present section we provide preliminary numeri- 
cal results: gap structure and density of states, obtained 
from the above theoretical formulation. 

The gap structure calculated by the procedure in § 2.1 
is shown in Fig. 2. Since we assume the orbital sym- 
metry k x ± iky is realized, the gap function possesses 
the in-plane fourfold symmetry (Fig. 2). We have ob- 
tained a gap structure which possesses strong in-plane 
anisotropy and band dependence. On the 7 band, we have 
a node-like structure to the directions [100] and [010]. 
This is because of the momentum-space periodicity and 
the odd-parity symmetry of superconducting gap func- 
tion, as pointed out by Miyake and Narikiyo. 18 On the 
small-gap bands a and 0, nodal structures are obtained 
on the diagonal directions [110] and [110]. This is because 
the p-wave attraction is weakened around the diagonal 
points by incommensurate antiferromagnetic fluctuation. 
This incommensurate antiferromagnetic fluctuation is at- 
tributable to the nesting between the Fermi surfaces a 
and 0. 45-48 Similar nodal structure on the a and bands 
was obtained even by using a two-band model for the a 
and bands. 49 

The density of states (DOS) is calculated by the for- 
mula 



cosh {yt) 



where I a (kF, z) is given by eq. (54). 



(64) 



1. 



PW) = --Im 

7T N 



(65) 



The density of states calculated for T = is shown 
in Fig. 3. The main part of the total DOS is taken by 
the main branch 7. The calculated partial DOS is 13.5 
%, 29.1 % and 57.4 % of the total DOS for the a, 
and 7 bands, respectively. This percentage is quantita- 
tively consistent with the expectation from the de Haas- 
van Alphen measurement. 38 The main coherence peaks 
around u = ±0.005 arc attributable to the DOS struc- 
ture on the 7 band. The low energy part of the DOS near 
the Fermi level is dominated by the small-gap bands a 
and 0. Our numerical result predicts that there are some 
fine structures between the main coherence peaks, which 
will originate from the coherence peaks on the a and 
bands. Actually, any spectroscopic measurements can- 
not detect such fine peak structures, due to their insuf- 
ficient resolution. However, some inflection points might 
be observed in p(u>) as a function of energy u by spectro- 
scopies. Such fine structures, if they are indeed observed 
experimentally, could be considered as the nature of the 
orbital-dependent superconductivity. 17 

3.2 Analysis of ultrasound attenuation rate 

Ultrasound attenuation rate is calculated by using 
the formula (53). The electron-phonon coupling matrix 
Afc ; q ; «' is calculated by using eqs. (43)- (46). The con- 
stants ijfj's are determined by fitting to the experimental 
results of ref. 14. In the present study we take gi = 0.192, 
.92 = 0.0096, g 3 = 0.0672, g 4 = 0.048 and g 5 = 0.03072 
for uj (q) = 1.0. The numerical results of ultrasound at- 
tenuation rate for various propagation and polarization 
directions are shown in Fig. 4. The remarkably strong 
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(a) 




Fig. 2. Calculated superconducting gap structure, (a) The gap 
magnitude on each of the three circular(cylindrical) Fermi sur- 
faces, a, /3 and 7 is depicted as a function of the azimuthal angle 
<j> around the c-axis (the Ru-O bonding direction corresponds to 
4> = 0). The unit of energy on the vertical axis is about 660(K). 
(b) The three Fermi surface sheets and the gap function. The 
Fermi surfaces are depicted by the thick solid lines. The depen- 
dence of the gap magnitude on in-plane direction is expressed by 
the distance from the Fermi circle along the direction. 

anisotropy is in semiquantitative agreement with the ex- 
perimental results (to be compared with Fig. 2 of ref. 14). 
The attenuation rate of the sound mode L100 is about 
one thousand times larger than that of the mode T100. 
As pointed out by Walker and collaborators, 34 this strong 
anisotropy is attributable to the anisotropy of electron- 
phonon interaction. 

In Fig. 5, the numerical data are compared with the ex- 
perimental data. In addition, the contributions from each 
band are separately presented there (note that the quan- 
tity is separable into the contributions from each band, 
since the formula (53) contains the summation with re- 
spect to band index a). 

The numerical results are well fitted to the experimen- 
tal results in the overall temperature region. For L100 
mode we can see that the ultrasound attenuation is dom- 
inated by that on the 7 band. In this case the contribu- 
tions from the a and (3 bands are almost negligible. On 
the other hand, for T100 mode, the sound attenuation on 
the 7 band is ineffective, and the contribution from the 
7 band is comparable to those from the a and j3 bands. 
The attenuation at low temperatures is dominated by 
the passive bands, particularly the a band. 
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Fig. 3. Calculated density of states in the superconducting state. 
The vertical axis represents the density of states (arb. units) 
and the horizontal axis represents the energy u> (the Fermi level 
corresponds to u> = 0). The thick solid line denotes the total 
density of states, and the thin solid lines denote the contributions 
from each band. The unit of energy on the horizontal axis is 
about 660(K). 

3.3 Analysis of thermal conductivity 

The thermal conductivity along [100] direction (the 
directions of the thermal current and the temperature 
gradient are both along the a-axis) is calculated using 
the formula (64). After the fitting procedure for specific 
heat, 23 there remains no fitting parameter, in contrast 
to the case of sound attenuation rate. The numerical re- 
sults are shown in Fig. 6. There we can show also the 
contributions from each band separately (note that the 
quantity is separable into contributions from each band, 
since the formula (64) contains the summation with re- 
spect to band index a). 

It is very interesting that actually the passive bands 
a and (3 contribute significantly to the thermal trans- 
port. Their contributions are comparable to that from 
the main branch 7. This situation is in contrast to other 
physical quantities in Sr 2 Ru04. For example, the main 
contribution to the specific heat is from the main branch 
7, 23 since the 7 branch takes the main part of the to- 
tal density of states. The reason why the main branch 7 
is not dominant in the thermal transport is as follows. 
The superconducting gap on the main branch 7 has a 
node-like structure near the zone boundaries (±7r, 0) and 
(0, ±7r). However, the Fermi velocity v kF around these 
points is quite small, and therefore the thermally excited 
quasi-particles around there can not play an essential 
role for the transport (note that the formula (64) con- 
tains multiplication of the Fermi velocity). This situa- 
tion demonstrates that the thermal transport property 
is inappropriate for detecting the gap anisotropy on the 
7 branch. 

4. Discussion 

In this section, we suggest some remarks from the 
present calculation. 

As we have seen in § 3.1, the gap structure is highly 
anisotropic. Although we have assumed naturally that 
the chiral state A(fc) ~ k x ±ik y is realized, the gap func- 
tion possesses sharp depressions on the (3 Fermi surface 
(Fig. 2). This anisotropy is sufficiently large for explain- 
ing the power-law temperature dependence of various 
physical quantities. We obtain a good fitting for spc- 
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Fig. 5. Calculated ultrasound attenuation rates (thick solid lines) 
are compared with the experimental data (o) read from ref. 14. 

(a) For the longitudinal wave propagating along [100] direction. 

(b) For the transverse wave propagating along [100] direction. 
The vertical and horizontal axes represent a(T) and tempera- 
ture T normalized by the value at T = T c , respectively. The 
contributions from each band are separately shown by the thin 
solid lines. 



Fig. 4. Calculated ultrasound attenuation rate for various sound 
propagation and polarization directions: (a) the longitudinal 
mode propagating along [100] direction and the transverse mode 
along [110] direction (i.e., L100: e//q, q//[100], and T110: e ± q, 
q//[110]), (b) the longitudinal mode along [110] direction (i.e., 
L110: e//q, q//[110]), and (c) the transverse mode along [100] 
direction (i.e., T100: e _L q, q//[100]). This figure is to be com- 
pared with Fig. 2 in ref. 14. 

cific heat, using the present gap structure, as we have 
shown in the previous work. 23 We should note that this 
is a new mechanism of creating nodal gap structure. 
The incommensurate antiferromagnctic fluctuation ac- 
tually causes the strong gap anisotropy on the a and 
(3 bands. Our microscopic theory suggests the follow- 
ing rule generally: singlet [triplet] superconducting gap 
function should change [not change] its sign between the 
Fermi surface portions which are bridged by antiferro- 
magnctic nesting vector, otherwise should take small val- 
ues (i.e., nodal structure) on these Fermi surface por- 
tions. Recently, Kontani has applied the same scenario 
to nickel-borocarbide superconductors, 50 which are con- 
sidered to be a s-wave superconductor having a highly 
anisotropic gap. He considered the effect of antiferromag- 
netic spin fluctuation on s-wave superconducting gap 
structure, and showed that the gap magnitude is sup- 
pressed in the Fermi surface regions which the antifcr- 
romagnetic nesting vector bridges. These results suggest 
that nodal gap structure does not always originate from 
the symmetry properties of superconducting order pa- 
rameter generally. So far, there have been a lot of phe- 
nomenological discussions. They have assumed that su- 
perconducting order parameter changes its sign on the 
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Fig. 6. Calculated thermal conductivity K[ioo] (T) (thick solid 
line) is compared with the experimental data (o) read from 
ref. 16. The vertical and horizontal axes represent «[ioo] {T)/T 
and temperature T normalized by the value at T = T c , respec- 
tively. The contributions from each band are separately shown 
by the thin solid lines. 

nodes, and have discussed the symmetry properties of 
order parameter on the basis of the temperature depen- 
dence of physical quantities. Such phenomenology has 
often been applied to uranium heavy fermion supercon- 
ductors. 51 However, as we have seen above, the validity 
of those phenomenological discussions is not so reliable 
in general. 

Now we turn our attention to the specific discussions 
on the gap structure of Sr 2 RuC>4. Recently Deguchi and 
Maeno measured specific heat under magnetic field, 52 ' 53 
in order to elucidate the detailed superconducting gap 
structure. They investigated the field-orientation depen- 
dence of specific heat. The anisotropy which is consid- 
ered to originate from the gap anisotropy on the 7 band 
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is observed to vanish at low temperatures. This might in- 
dicate the possibility that the small-gap bands (a and (3) 
possess nodal structures near the diagonal lines, and the 
in-plane anisotropy from the active 7 band is canceled 
by that from the passive a and (3 bands at low tem- 
peratures. 53 Recently Kusunose has analyzed the field- 
orientation dependence of H C2 and specific heat C. 54 
He showed that the gap structure with the intermediate 
magnitude of minima in [100] direction for the 7 band, 
and tiny minima of gaps in [110] directions for the a and 
(3 bands give consistent behaviors with experiments. Par- 
ticularly, he succeeded in explaining the anomalous tem- 
perature dependence of the in-plane anisotropy of H C2 , 54 
i.e., the sign change of SH C2 near T c . 55 These experimen- 
tal and theoretical results are consistent with our gap 
structure. Further investigations should be performed to- 
ward the thorough clarification of the gap structure of 
Sr 2 Ru0 4 . 

We have used the electron-phonon coupling constants 
(ji's as fitting parameters. It is very difficult to deter- 
mine these parameters microscopically without the fit- 
ting. In order to determine the parameters gt's without 
the fitting, it will be indispensable to know microscopi- 
cally how the atoms are rearranged in order to relax the 
local stress (or strain), because g^s will essentially de- 
pend on the way of the rearrangement. For example, the 
stress along [100] direction can be relaxed by reducing the 
Ru-O-Ru angle from 180 degrees without shrinking the 
Ru-0 bonds along [100] direction, as well as by shrinking 
Ru-0 bonds without changing the Ru-O-Ru angle, cji's 
will be different between these two ways of relaxation. 

We have found actually that the overall temperature 
dependence and anisotropy of the ultrasound attenuation 
rate depend significantly on the parameters g^s. For ex- 
ample, we have shown the results in which a(T)/a(T c ) 
for L100 is sma/Zerthan that for T100 (Compare Fig. 5(a) 
with Fig. 5(b)). However, we can actually show by using 
other parameter sets of g^s that a(T)/a(T c ) for L100 
could be larger than that for T100. Therefore it is ac- 
tually difficult to determine the superconducting gap 
anisotropy by the experimental data of ultrasound at- 
tenuation rate without any detailed information on the 
electron-phonon coupling matrix elements. 

We find still slight deviation between the theoretical 
and experimental results of ultrasound attenuation rate 
in Fig. 5, particularly at low temperatures in Fig. 5(a). 
We could consider various reasons for the deviation. One 
possible reason is that we could not reproduce the Fermi 
surface deformation precisely by using only five compo- 
nents of the electron-phonon coupling matrix elements. 
To obtain better agreement, we might require higher or- 
der harmonics of the matrix elements. 

Recently Contreras et al. have obtained a good fitting 
of the ultrasound attenuation rate in Sr 2 Ru04. 3 ° Note 
that their gap structure is quite different from ours. They 
claim that point nodes should be located on the 7 and 
(3 bands. 30 However, it might be difficult to reproduce 
the line-node-like behaviors observed experimentally in 
many quantities. In addition, since the jump of specific 
heat at T c will be larger than that in the line-node case, 
it is unclear whether or not to obtain consistency with 



the experiment of specific heat. From our microscopic 
point of view, we expect that, since the pairing interac- 
tion could possess only negligible momentum dependence 
along the c-axis due to the strong two-dimensionality, 
point nodes as well as horizontal line nodes are almost 
impossible to realize in Sr2Ru04. 

The ficld-oricntation dependence of thermal conduc- 
tivity has been utilized to investigate the anisotropy of 
superconducting gap. 15,16,56 Izawa et al. observed actu- 
ally negligible in-plane fourfold-symmetry component in 
the field-orientation dependence of the thermal conduc- 
tivity, and concluded that line node should run around 
the cylindrical Fermi surfaces (horizontal line node). 15 As 
we have shown in the present work, thermal conductiv- 
ity is ineffective for detecting the gap anisotropy on the 
7 band. This will be the reason why the gap anisotropy 
of the 7 band can be observed in field-oriented specific 
heat, 52 ' 53 but not in field-oriented thermal conductiv- 
ity. 15 ' 56 Therefore one might expect that thermal con- 
ductivity is appropriate for detecting the gap anisotropy 
on the passive bands a and (3. However the fourfold- 
symmetry component of the field-orientated thermal con- 
ductivity may barely be detectable for the a and (3 bands, 
as Tanaka et al. demonstrated. 27 We consider that the 
field-oriented thermal conductivity does not crucially ex- 
clude the possibility of line nodes along the c-axis. 

In any case, we could conclude that line-node-like 
structure should exist on the a or (3 band, by combin- 
ing the following two points: (1) the T-linear behavior of 
k/T indicates that line-node- like structure should exist 
somewhere on the Fermi surface and (2) the thermal ex- 
citation on the 7 band does not effectively contribute to 
the thermal transport. However, it would still be difficult 
to obtain some insights into the nodal position only from 
the present analysis of thermal conductivity. 

Recently Yanase et al. proposed another type of gap 
structure. 57 Although the essence of their pairing mech- 
anism is the same as ours, 58 the gap structure is a little 
different from ours: the nodal position on the (3 band de- 
viates from the diagonals. It is still controversial whether 
such a deviation is realized or not. 

For determining the nodal positions thoroughly, it 
would be required to measure the gap magnitude by some 
momentum-resolving experiment on the Fermi surface. 
Further experimental probes crucially determining the 
nodal positions are desired to be developed. 

5. Conclusion 

In the present article, we have discussed the transport 
properties in the superconducting state of Sr 2 Ru04, by 
analyzing the sound attenuation rate and thermal con- 
ductivity. We have shown that the gap structure given 
in Fig. 2 is consistent with the experimental data. 

As we have seen in Fig. 2, we obtain a nodal gap struc- 
ture with a large in-plane anisotropy, particularly on the 
passive bands a and (3, even if we assume the chiral state 
A(fc) ~ k x ± iky. We should note this is a new mecha- 
nism of creating nodal gap structure: antifcrromagnctic 
fluctuation deforms the gap structure, and creates the 
nodal gap structure. We do not consider that such a 
nodal structure due to this mechanism occurs exception- 
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ally only in Sr 2 Ru04. We consider that a nodal struc- 
ture could occur due to this mechanism also in other su- 
perconductors, particularly likely in strongly correlated 
electron systems. Therefore we should notice generally: 
actually it is not reliable to determine pairing symmetry 
by power-law temperature dependence of physical quan- 
tities. 

Through the present study we would like to stress that 
too much simplified models (e.g., an isotropic Fermi sur- 
face, only a single band or a gap function approximated 
by only a few harmonic functions.) are insufficient in gen- 
eral to analyze the anisotropy and the temperature de- 
pendence of physical quantities. In addition, we should 
note in general that some quantities can succeed in de- 
tecting the anisotropy of superconducting gap, but oth- 
ers can not. An example is the thermal conductivity in 
Sr 2 Ru04, which would fail to detect the in- plane gap 
anisotropy on the 7 band, as we have discussed in § 4. 

In any case, we consider that it is very difficult to ex- 
plain the experimental data without assuming the exis- 
tence of line nodes on the passive small-gap bands a and 
p. 

The present work is one of a few examples in which 
physical quantities are analyzed with use of realistic 
many-band electronic structure and superconducting gap 
structure obtained by microscopic calculation. In conclu- 
sion, we hope that such detailed analyses are performed 
also for many other anisotropic superconductors. 
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Appendix A: Electron-Phonon Coupling Ma- 
trix 

We consider the electron-phonon interaction, since the 
change of electronic structure due to lattice distortion 
is essential for the ultrasound attenuation. The detailed 
information about the anisotropy of electron-phonon in- 
teraction is actually important for the realistic analysis, 
as Walker et al. pointed out. 34 

The non-interacting part of the Hamiltonian Ho is 
written in the form, 

H = ^2 t ee'[ r i- r i']4e* c i't'<r- (A' 1 ) 

ii',M',<T 

tu'[ r i — r i'\ is the transfer matrix element between the 
Wannier atomic orbitals I at i-th site and £' at i'-th site 
. Now the lattice deformation is introduced: 

H = ^2 + u i) - ( r *' + u v)\ c \ic c v^<y (A-2) 

Ui is the displacement of the i-th Ru atom. The hopping 
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integrals are expanded in uf. 

t U '[(ri+Ui) - (tv +«.')] 
= t U '[ri - tv] +g u >(ri - tv) • (m - u^) + (A-3) 

where g u ,(ri — r^) — dta>[a]/ da\ a = ri - r ., The Hamilto- 
nian of the electron-phonon interaction is given by 

H cp = 9w( r i - r i') ■ ( u i - uv)4io c VL'°- (A-4) 

Then we transform the expression to the momentum rep- 
resentation, and apply the second quantization to the 
lattice vibration, 

ffc P = J2 (9ee'(R)-e)J2( 2NM ^(q)r 1/2 

R,U',a q 

x ^ ife '*(l - e- { o R ){b q + bl q )c{ +qlcj c kfJ „. (A-5) 
k 

Thus the Hamiltonian takes the form of eq. (42), using 
the electron-phonon coupling matrix 

A fc ,,,«, = (2Af W „(<j))- 1 / 2 ^e- iM 

R 

x(g u ,(R) ■ e)(l ~ e~ iq - R ). (A-6) 

M is the ionic mass, u q is the phonon frequency, 
and b q (b q ) is the annihilation(creation) operator of the 
phonon with momentum q and polarization e. The unit 
vector e points to the phonon polarization direction. In 
the actual experiment, the phonon wavenumber is much 
larger than the lattice constant (q <C 1), while the sum- 
mation in R of eq. (A-5) is performed at most up to the 
second nearest sites. Therefore, H ep is approximated as 
follows by expanding the factor e~ lq R in the power of q 
and neglecting higher orders: 

Ak, q .u' = i(2Mcj (q))- 1/2 ]T e- ik - R (g u ,(R) ■ e)(q ■ R). 

R 

(A-7) 

Here we consider the specific expression of the function 
g U i{R). Since we expect naturally that the electronic 
transfers perpendicular to the direction of the deforma- 
tion are not significantly affected by the deformation, we 
may assume that g u ,{R) is parallel to the vector R: 

g u ,{R)=g u ,{R)R. (A-8) 

R is the unit vector along R, i.e., R — R/\R\, 
and ga'(R) is an even function of R [i.e., gw(R) = 
gw(-R)]- The coefficients gw{R) are characterized by 
five constants, since we have taken the five hopping in- 
tegrals (ti, ...,t$) for describing the electronic structure: 
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Fig. B-l. The diagrammatic representation of the lowest-order 
electron self-energy correction due to the lattice distortion r]<j> q ■ 
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(A-16) 
(A-17) 
35, (A-18) 



and the other elements of gw(R) are zero, where x(y) 
is the unit lattice vector along the x-(y-)axis. We obtain 
the momentum-dependent electron-phonon coupling ma- 
trix elements in eqs. (43)-(46), by substituting eqs. (A-9)- 
(A-18) into eq. (A-7) and using the coefficients g^s re- 



lated to <7j's as gi — 2^/gi, g 2 = 2 3 /' 



732, 33 



2733, 



54 = 2 7ff4 , 35 = 2 3 / 2 735 (7 - {^o(q)/(2M W 2)}i/2 ; and 
v q is the sound velocity). 

Appendix B: Charge Neutrality Condition un- 
der Lattice Deformation 

In this Appendix, we consider the local charge neutral- 
ity condition, which should be retained under the lattice 
distortion due to the ultrasound propagation. We intro- 
duce a lattice distortion described by rj(j) q [For the def- 
inition of <j) q , see eq. (42)], whose spatial periodicity is 
characterized by the wavenumbcr q. r\ is related to the 
magnitude of the displacement of the atoms. This dis- 
tortion plays a role of periodic potential acting on the 
electron system. The electron self-energy will be modi- 
fied by 



H u ,{k;q)=r 1 D{q)K 



k,q,W 



(B-l) 



Now we focus on a local small volume, which is enough 
smaller than the characteristic length of lattice distor- 
tion (~ 27r/|q|) but still contains macroscopic number of 
electrons. In this volume, the electron system is regarded 
as spatially uniform. As far as we discuss the properties 
of the electrons contained in this small volume, we can 
use the Green's function Gw{k) corrected by the self- 
energy T,££i(k;q). The local change of electron number is 
obtained by 

T 
l N 



Sn 



2^^[G«(fc)-G^(fc)] 



a.k 



within the precision to the lowest order in rj, and we have 
used similar relations to eqs. (23) and (24). The charge 
neutrality condition states that the charge inhomogcne- 



ity caused by the lattice distortion rj<j) q should be can- 
celed by introducing a local change of electron chemical 
potential: 41 

1 5> ? G%(fc) \^ ul fi {k)U t2a {k)^ lt2 {k ]q ) 



a.k 



where (d^G) a (k) is defined as 

(d ( G) a (k)^Tj2 



We obtain from eq. (B-3) 



-A/i 



dd 0) (fc) 



0, (B-3) 



(B-4) 



A/t = 



E a , k (d 6 G) a (k) 

(B-5) 

Here let us note that the shift of chemical potential 
A/x can be regarded as a renormalization of the electron- 
phonon coupling matrix elements Ak,q,ll'- The chemi- 
cal potential shift A/i appears in the equations always 
together with the self-energy in the form, T,w(k;q) — 
Afidw. If we define the mean value of \k,q,ee' by eq. (56), 
then 



Y,u>{k;q) - A/z5«< = rjD(q)(A k . 



q,U' 



0- (B-6) 



Thus the local charge neutrality condition is satisfied by 
using the renormalized value of the electron-phonon cou- 
pling, Ak, q ,U' - A q Su>, instead of Ak, q ,w- 

In the normal state, using eqs. (11) and (B-4), we ob- 
tain 



(d s G) a (k) = 



<9np(e) 



de 



(B-7) 



=l«(fe) 



riF(e) is the Fermi function. In the superconducting state, 
we obtain eq. (57) by using eq. (B-4) and Gi 0) (k) of the 
form 

iw„ +£ a (fc) 



GW(fc) 



Uk) +\^a(k) 



(B-8) 



Appendix C: About the Strong Gap 
Anisotropy on the a and (3 
Bands 

In this Appendix, we discuss the origin of the strong 
gap anisotropy on the a and bands by extracting the 
pairing interaction V ar7lCr2ia ' a3Cr4 (k, k') in cq. (12). We fo- 
cus on the values of Vg a(Jy saa{k,k') and Vg aayaa(T (k, k') 
on the Fermi surface. We characterize position on the 
Fermi surface by the in-plane azimuthal angle <fi with re- 
spect to the Ru-0 bonding direction. Hereafter we refer 
to the values of Vg aa jj aa (k, k') and Vg acr . acrcr (k, k') on the 
Fermi surface as Vg : g(<j)g; (j)' g ) and Vg^^g;^), respec- 
tively (The Matsubara frequencies u> n and w' n are fixed 
to ttT). Now we fix the initial state of Cooper pair on 
the Fermi surface by <j>g = \ as pointed in Fig. C-l(a). 
Vp,p{<pp = 7r/4; <f>' 3 ) and Vp t0t (<j>p = tt/4; <j>' a ) are depicted 
as a function of the angles 4>' g and <f>' a , respectively, in 
Fig. C-l(b). We find that there are characteristic local 
minimum in Vp t p{(f>p = it / 4; (f>' g ) around (j)' g — 57r/4 and 



12 J. Phys. Soc. Jpn. 



Full Paper 



Takuji Nomura 




<t»'tf=7l/4 




















































i|>'p=5jc/4 



Fig. C-l. (a) The three Fermi surfaces, a, (3 and 7 are depicted. 
Position on the Fermi surfaces /3[a] is characterized by the in- 
planc azimuthal angle <j> around the points r(0, 0)[X(-7r, w)] with 
respect to the i-axis. (b) Calculated effective pair-scattering am- 
plitudes Vp t p(<j>p = ir/4;rf>'p) and Vp ta ((fif3 = n / 4; <j>' a ) are shown 
as a function of the azimuthal angles <j>L and <j>' a , respectively. 



local maximum in Vf3.a(<t>/3 — tt/4; 4>' a ) around <f>' a — n/4. 
The local minimum of Vp t /3(<j>p — 7r/4; (/>'*) around <p'g — 
5ir/4 favors the same sign of order parameter between the 
points (j>p — 7r/4 and (f>'p — 57r/4. The local maximum of 
Vp,a(4>i3 ~ 7r /4; 4>' a ) around <f>' a — ir/4 favors the different 
sign of order parameter between the points tfip — ir/4 
and <j>' a — 7r/4. These momentum dependences deform 
the gap structure on the a and (3 bands significantly, and 
result in the strong anisotropy of gap structure. Those 
characteristic momentum dependences of the pairing in- 
teraction basically originate from those of the fluctuation 
due to the nesting between the a and f3 Fermi surfaces, 
although it is somewhat modulated by the diagonaliza- 
tion matrix U(k) in eq. (10). 
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